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The power of linearity and superposition I 
The time domain: Convolution and the Impulse Response 

Convolution can be used to calculate the response of a system to arbitrary inputs by using the impulse response of a 
system. 

As an example consider the system described by the differential equation: 

y(t)+y(t) +y (t) = f(t) 

where y(t) is the output, and f(t) is the input. The impulse response of this system, h(t), can be shown to be equal to 
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The impulse response is shown below. 
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You have learned techniques for finding y(t) if f(t) is a well defined input such as a step, impulse or sinusoid. 
Convolution allows you to determine the response to more complex inputs like the one shown below. 
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There are several ways to understand how convolution works. First convolution will be developed in an approximate 
form as the sum of impulse responses. This presentation is useful for an intuitive understanding of the convolution 
process. After the approximate form is developed, the exact analytic form of convolution is given. You can also link to 
an example of the convolution integral in action , and to a less physical (more mathematical) derivation . 



The convolution as a sum of impulse responses- 

(the Matlab script, Convolution.m , was used to create all of the graphs in this section). 



To understand how convolution works, we represent the continous function shown above by a discrete 
function, as shown below, where we take a sample of the input every 0.8 seconds. 




Input vs Time and discrete approximation 
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The approximation can be taken a step further by replacing each rectangular block by an impulse as shown 
below. The area of each impulse is the same as the area of the corresponding rectangular block. 



Input vs Time and approximation with impulses 
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The superposition theorem states that the response of the system to the string of impulses is just the sum of 
the response to the individual impulses. The response of the system to the individual impulses is shown 
below. 
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Lets look carefully at this graph and note the following: 

• The response of the system to a unit impulse, 8(t) , is h(t). 

• The response of the system to a unit impulse delayed by time i ■ AT , 8(t - i ■ AT) , is h(t - i ■ AT) . 

• The response of the system to an impulse of area A delayed by time i ■ AT , is A ■ h(t - i ■ AT) . 

• The height of the function at a time t = i ■ AT is f (i ■ AT) . 

• The area of the impulse at t - i ■ AT is f (i ■ AT) ■ AT . 

• The delayed and shifted impulse response is given by f (i ■ AT) ■ AT ■ h(t - i ■ AT) 

• Therefore the response of the system to a series of delayed impulses, as shown above is 

y(t)=jr(f(i-AT)AT)h(t-i-AT) 

In this equation 

o y(t) is the output 

o i ■ AT is the time delay of each impulse 
o (f (i ■ AT) - AT) i s the area of the I th impulse 

o if you take the limit as AT^> 0 ? the summation yields the convolution integral (with 
i' AT = r, AT = At) 

y(t)=J o f(T)h(t-T)dT 

This is the Convolution Theorem, 
o The integral is often presented with limits of positive and negative infinity: 

.y(t)=£f(T)h(t-T)dt 

For our purposes the two integrals are equivalent because f(t)=0 for x<0 
o The form of the integral we will most often use is 

y(t)=J o t f(T)h(t-T)dT 
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This form is also equivalent because h(t-x)=0 for t-t<0 (or h(t-t)=0 for t>t). This makes the 
limits of integration effectively 0 and t. 
o The arguments in the integral can also be switched to give two equivalent forms of the 
convolution integral 



3<t)-£f(tW-*)*=0<*)f(t-t)*; 



This equation merely states that the output is equal to the sum of the responses from the individual impulses. 
Another (more mathematical) derivation of the convolution integral is given here . 

The summed response is shown in the graph below, along with the individual responses. 
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Obviously if you decrease the accuracy of the approximation should improve. The graph below shows the 
impulse responses if we let AT=0.3 seconds (instead of 0.8 seconds as shown above). 
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The graph below shows our approximate response, for AT=0.8 and AT=0.3, along with the exact response 
(the "exact" response is actually a numerical approximation, as calculated by Matlab). 
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The graph below shows the "exact" response along with the response calculated using Matlab 

"conv" (convolve) function. Type "help conv" at the Matlab prompt to see how it works. Examine the script, 

convolution.!^ which generated all these graphs for an example of how to use the function. 




Hopefully this background gives you some physical insight into how convolution can be used to find the 
response of the system to an arbitrary input. 



An analytic description of convolution 

To describe the analytic form of convolution consider the convolution integral 
y(t) = rf(T)h(t-T)di 
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or equivalently 



y(t)=J t f(T)h(t-T)dt 

if h(t)=0, t<0 and f(t)=0, t<0. 

Consider the system discussed previously with the impulse response described by: 
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as shown below 
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To find the response of the system, y(t), to an input, f(t), at t=l sec we need to evaluate the integral 



y(t)=p(t)h(t-T)dT 



with t=l. To visualize this process consider the product f(x)h(t-x). A plot of f(x) yields the input. A plot of h 
(t-x) yields a time-reversed version of the impulse response with a time delay of t (1 second in this case). 
These are shown below along with the product of the two functions (solid black line). The function h(t-x) is 
plotted in magenta and doesn't look much like the function above -- but note the different vertical scales. 
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The next graph shows the result of the integration for all time, with a black dot at t=l . 
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As t is increased further the function h(t-x) shifts to the right. The multiplication and integration can be 
evaluated for each value of t to yield the system response y(t). The graphs below are similar to the ones 
above, except that t=4. 
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To see an animation of this shifting, multiplying and integrating process, click here. 

Using the convolution integral it is possible to calculate the output, y(t), of any linear system given only the 
input, f(t), and the impulse response, h(t). However, this integration is often difficult, so we won't often do it 
explicitly. Later you will learn a technique that vastly simplifies the convolution process. 



Some more examples 

To see some more animations of convolutions, click here . 
Another demonstration 

Another excellent web-based demonstration of convolution is The Joy of Convolution at Johns Hopkins. 
An Example of the Convolution Integral in Action 

This page has given a description of the convolution process, but has not actually gone through the 
mathematical procedures needed to analytically evaluate the convolution integral. This process is described 
on another page. Click here to go to that page . 



The MatLab code that generated the animations is available here . 



